Topologically induced vortex magnetance of a quantum device 
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Abstract 

We consider resonant vortices around nodal points of the wavefunction in 
electron transport through a mesoscopic device. With a suitable choice of the 
device geometry, the dominating role is played by single vortices of a preferred 
orientation. To characterize strength of the resulting magnetic moment we 
have introduced a magnetance, the quantity defined in analogy with the device 
conductance. Its basic properties and possible experimental detection are 
discussed. 
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The topological structure of quantum mechanical wavefunctions is responsible for many 
observable phenomena. One of its prominent consequences is existence of vortices. They 
have already been observed in macroscopic quantum systems, rotating superfluid helium and 
superconductors. In general, their existence is connected with the time reversal symmetry 
breaking. The first example of the such topological effect was already given in the very early 
days of the quantum theory. It concerns the probability current density of a single electron 
defined conventionally by 



which may be nontrivial once the wavefunction tp is complex. In systems with open ge- 
ometry a nonzero probability current refers to an electron transport through the system, 
i. e. from the source to a drain. DiracS pointed out that the quantum probability current 
may exhibit a vortex structure around nodal points (zeros) of the corresponding wavefunc- 
tion. The effect was later discussed by Hirschfelderi. Recently it has been demonstrated 
numerically that the current in topologically nontrivial devices exhibits pronounced vortices 
the form and magnitude of which changes quickly with the energ yll and applied magnetic 
fieldi. The experimental evidence of their existence is of the particular interest. 

The direct consequence of a vortex is a non-zero magnetic moment, which has to depend 
on the applied current in a similar way as the device conductance. The natural physical 
quantity representing a vortex structure strength is thus the device magnetance given by 
the ratio of the magnetic moment and applied voltage drop. To analyse the effect of vortices 
to the electron transport and magnetance we will limit our consideration to devices with 
two-dimensional gas of spinless electrons. Main attention will be payed to the simple device 
formed by a piece of the one-dimensional wire with tangentially attached circular cavity (a 
quantum dot) as sketched in Fig. 1. It may be expected that there appears a dominating 
circulating current in analogy with the classical gas flow. This structure is thus the natural 
candidate for an experimental attempt to observe electron vortices. 

The vortices are closely related with the topological structure of the phase. To be more 
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explicit, let us write the wavefunction as 

1> = y/pe iS , (2) 

where p is the probability density of the particle and S is the corresponding phase. The 
latter is defined modulo 2tt and assumes conventionally values in the interval [0, 2n) . In 
the two-dimensional case it is possible to identify the configuration space with a region in 
the complex plane and to treat a possible multivaluedness of S in terms of its analytical 
structure, especially cuts and branching points corresponding to different Riemannian sheets 
of this function. 

Existence of a branching point on the phase Riemannian surface implies appearance of 
a vortex in the corresponding quantum probability flow!. Inserting (0) into ([!]) we get 

j(r) = A p(f) v5(f). (3) 

The vector v = (K/m*)vS can be regarded as a velocity of the corresponding probability 
flow. For a closed curve T the vorticity of v along T is thus 

/ v dl = —5S = 2vr— m ; m = 0, ±1, ±2, . . . , (4) 
Jv m* m* 

where 5S is the phase change when winding once around the curve. Since the wave function 
(@) must be single- valued, the difference 5S can only equal to a multiple of 2n . If the phase 
S has no singularities inside T, the contour can be shrinked into a point in which case the 
vorticity is zero. If, on the other hand, T encircles a nodal point of if), the phase is ambigu- 
ous at this point and the integer m may be nonzero. In this situation the corresponding 
probability current exhibits a vortex centered at the nodal point. 

In quantum devices three basic topological situations can occur (see Fig. 1) : 
(a) A phase cut starts and ends at the boundary of the system. This is typical for integrable 
systems. As a simple example, consider the transport through a straight quantum waveguide 
of the width w parallel to the x-axis. If the incident wave is in the n-th transverse mode, 
the wavefunction has the form 
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ipk,n{ x ,v) = \ —e ikx sin — , (5) 

V 7TW W 

where n is the mode number and h 2 ( k 2 + 7r 2 n 2 /w 2 )/2m* is the particle energy. In this case 
the phase S = kx is monotonically increasing in the longitudinal direction and its cuts are 
located at the segments x = 2irj /k , j = . . . , —1, 0, 1, ... , parallel to the x-axis. Their 
endpoints lie at the boundary, y — 0, w . They are physically irrelevant being not branching 
points. 

(b) The cut starts at the boundary and ends at some nodal point inside the system. The 
internal endpoint is a branching point which is related to a single vortex. 

(c) The cut connects two nodal points of the wavefunction. Both endpoints are branching 
points which corresponds to a pair of vortices rotating in opposite directions. 

To define the vortex magnetance we will follow the scattering approach!, generally ac- 
cepted in the transport theory of quantum devices. To establish conductance of a two 
terminal device, infinitely long ideal leads are placed between the device (scattering region) 
and electron reservoirs (source and drain) allowing an explicit asymptotic form of scattering 
wave functions. The obtained transmission coefficient T(E) is used to relate the applied 
current J and the chemical potential difference between reservoirs, Afi. In the limiting case 
of vanishing Afj, and at zero temperature the electron transport is determined by properties 
of electrons at the Fermi energy Ep and we have 

J = ^ T(Ep) Afi . (6) 

The applied current is responsible for a non-zero momentum of the system. It can be 
divided into two parts: momentum of the center-of-mass and the momentum due to electron 
motion relative to the center-of-mass. The later one originates in circulating currents giving 

— # 

rise to a magnetic moment M = (0, 0, M) representing their strength. Momentum of the 
center-of-mass is controlled by the current density within the ideal leads only. Note that 
in the considered limit of infinitely long ideal leads the mass center momentum cannot be 
affected by a device of finite dimensions. 
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In the simplest case of one-dimensional ideal leads having a common axis the momentum 
of the center-of-mass can easily be evaluated. It is controlled by the product T(E F )j (r) 
with jo(r) being the current density within the system without the scattering region (the 
device is replaced by the lead). In this case we get the vortex magnetance in the following 
form: 

M(E F ) eg(E F ) r n 2 

Aj2 = 2c J (pin -T[E F )j (r)) xrd r, (7) 

where g(E) denotes the density of device states. The magnetance defined by this way 
characterizes the vortex structure of the studied device and it is invariant with respect to 
coordinate system translations. 

The expression, Eq.([7|), is applicable to the device sketched in Fig. 1. To stress the effect 
of the device geometry the flat potential is assumed within the device area demarcated by 
hard walls. The relevant wave functions are thus eigenfunctions of the Schrodinger equation 
for free electrons with zero values at the boundaries of the device and ideal leads. 

To estimate the magnetance of the considered device, let us first summarize wavefunc- 
tion properties of the separated circular cavity. In polar coordinates eigenfunctions are 
determined by zeros of Bessel functions J m at the cavity radius R and we have 

tf m (r,0) = ct, m J m (h m r)e ±imd , (8) 

where k^ m = x^ m /R with x^ m being the Z-th zero of the Bessel funct ion Jmi and ci m 
is the normalization factor. The quantum number m (m = 0, 1, 2, . . . ) represents the 
angular momentum ±(h/m*)m which is closely related to the vorticity defined by Eq.([|). 
The eigenstates of spinless electrons are two-fold degenerated with respect of the sign of the 
momentum. 

Only the azimuthal component of the probability current density for given eigenenergy 
Ei >m might be nonzero and it is given as follows 

Tl 771 

je{r) = A hm — —cl m Jl l {ki, m r) , (9) 
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where A^ m is given by the difference of weight factors, A^ m = |a|^ | 2 — K ^ | 2 , representing 
amplitudes af^ and a\ J of eigenfunctions with positive and negative values of the orbital 
momentum, respectively. In the equilibrium case the problem has time reversal symmetry 
requiring equality of weight factors and the total momentum vanishes. 

Attaching a wire the wavefunctions ip(r,6) given by Eq.(^) become modified by a cou- 
pling with plane waves. Assuming a small window between wire and cavity, much less than 
R, wave functions ip(r, 8) will only be slightly modified. The main effect will be a level broad- 
ening represented by the spectral density gi jUl (E) with sharp maxima at energies E = E^ m . 
Wave functions representing states with the same direction of the probability flow in the 
window region will be much easily matched together than those with opposite flow direc- 
tions. In the current currying regime the amplitude of cavity wave functions matched with 
incident plane wave will be thus enhanced giving rise to circulating currents. Corresponding 
magnetance may be estimated by choosing a proper value of the A^ m (E) entering Eq.(^). 
Inserting the current density j$(r), Eq.(|9|), into expression for the magnetance, Eq.(^), we 
get 



M(E F ) 



H* B J2 A i,m(E F )mgi :rn (E F ) , (10) 



where fi* B = eh/2m*c is the effective Bohr magneton. Due to the dominating role of the 
spectral density gi jm (E) a resonance character of the magnetance is expected. Since gi tTn (E) 
is proportional to m* the magnetance is of the purely topological origin. 

As an example, the energy dependence of the magnetance and the transmission coefficient 
close to the resonance for / = 1 and m = 3 is shown in Fig. 2. Corresponding current density 
distribution for Ep indicated by arrow is presented in Fig. 1. Four phase-cut endpoints 
represent vortex centers. Three of them are located close to the cavity center giving rise 
to pronounced circulating current. For fixed window region the broadening of cavity levels 
becomes weaker for larger radius R giving rise to more pronounced resonances as seen in 
Fig. 3a. For R > 1.8 w more complicated resonance structure appears due to states with 
different quantum numbers I. 



For a quite well separated resonance the magnetance peak height may be estimated by 
equality of the current originated in the cavity eigenf unction, Eq.(^j), and that of the single 
wire, Eq.(|5]), in the common region, i.e. in the area of overlapping circular and strip regions. 
This condition gives: 

hg(E r ) f R 2 2 dr 



A hm (E r ) —m / Cl J m (ki >m r) — 

m* JR-d r 



= - — / sm — dy , (11) 

h 7VW Jw-d w 

where d is the width of the common region, d = w + R — yo with y being the distance of 
the cavity center from the bottom wire edge. This estimation is in good agreement with the 
values obtained by direct evaluation of the Eq.(^), as shown in Fig. 3a. 

The above amplitude estimation, Eq.(0), assumes an ideal coupling allowing an easy 
transfer of the incident wave with transmission coefficient T approaching unity. Out of the 
ideal coupling a reflection will take place and T will decrease. However, far enough from 
the resonance the mixing of plane waves with eigenfunction of the separated cavity becomes 
less effective and transmission will be enhanced. This typical energy dependence is shown 
in Fig. 2b. 

At finite temperatures the magnetance is given by the following expression 
M(/i,T) 1 rdfa(E) 



, i M{E)dE , (12) 

Afi Afi J dp. 

where fo(E) denotes the Fermi-Dirac distribution function and M(E) is the zero tempera- 
ture magnetic moment given by Eq.(|7]). The resulting magnetance at a non-zero temperature 
is shown in Fig. 3b. For high enough temperatures resonances will fully be smeared out and 
the magnetance should approach values obtained by a semiclassical treatment. Since the av- 
eraged density of states g scales with the cavity area < g >= (m* /hh)Ep vri? 2 and m ~ R/w, 
the approximate scaling of the magnetance with R 3 is expected. 

The semiclassical description of the studied ballistic transport can be based on the billiard 
mo defi Electrons are injected uniformly over the wire width w and with the probability 
| cos cf) along the direction represented by the angle <fi with respect of the wire axis. As 
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expected, the resulting magnetance, shown in Fig. 3b, approaches mean values obtained by 
full quantum-mechanical treatment. 

The same qualitative features of the magnetance are expected for any cavity defined by 
a potential of the circular symmetry. We have found that resonances are very stable with 
respect of the radius modulation (AR(9)/ (R) <0.1) and potential fluctuations less than the 
energy level separation of the unperturbed cavity. Interaction between electrons, which was 
not taken into account, should not change the magnetance structure qualitatively. However, 
we expect enhanced magnetance due to the electron-electron friction. 

Impurity potentials within wire area have destructive effect to the magnetance. They 
suppress the transmission coefficient and consequently the current flow along window be- 
tween wire and cavity, which is responsible for the magnetance peak height as discussed 
above, Eq. (|TT|) . They also give rise to additional vortices of the uncontrollable orienta- 
tion within the wire region. The resulting richer magnetance structure can be much easily 
smeared out by the temperature. This destructive effects will be even more effective for 
wider, multimode, wires. 

In real systems the probability to have impurities within wire region might be lowered 
by the shortening of the wire length. Similar device, but with centrally attached wires rep- 
resented by point contacts, have already been realized^. The proper design could also allow 
to control cavity area by top gates. To observe resonances the energy level separation of 
the isolated cavity has to be larger than the thermal energy k^T . For the device fabricated 
from a GaAs-AlGaAs heterostructure it limits the cavity radius, -R[/im] < 0.2/yT[K]. For 
example, assuming R ~ 0.5/im resonances could be indicated by a fluxmeter allowing mea- 
surements with accuracy 10 2 fiB, which is in principle within experimental reach@. The 
contribution of non-circulating currents cannot be easily excluded by experimental setup 
as has been done in our theoretical treatment. Under the standard conditions of the fixed 
applied current it will be, however, responsible for a monotonic background only. 
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FIGURES 

FIG. 1. The model device geometry with phase cuts (solid lines) for the energy indicated in 
Fig. 2 by the arrow. The window between wire and the cavity equals to 5/4 of the wire width w 
and the cavity radius R = 5w/3. Corresponding current density is illustrated by arrows. 

FIG. 2. Energy dependence of the magnetance M/Afi (a) and transmission coefficient T (b) in 
the vicinity of the resonance for / = 1 and m = 3. Device geometry is the same as that in Fig. 1 
and E\ denotes the energy of the lowest transversal mode E\ = {ir 2 U 2 ) j '(2m* w 2 ) . 

FIG. 3. Magnetance M/A/j, at zero temperature (a) and at the temperature 0.05 Ei/ks (b) as 
function of the cavity radius for fixed window width bw/4. The dashed line represents estimation 
for resonance maxima, Eq.(ll). The dashed-dotted line is the result of the semiclassical treatment. 
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